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ABSTRACT 


ill 


Recently there has been considerable interest shown 
in Nuclear Reactors of newer desipis with sinplified 
features, These have sought to improve the safety aspect and also U 
the plant construction costs by sirrplifying features of the 
conventional reactors. The Natural Circulation Boiling Water 
Reactor is one such design which loc*;s very promising from the 
point of view of future power requirements. This thesis analyses 
the feasibility of such reactors satisfying one of the criteria of 
thermal hydraulic requirement, namely the Critical Heat Flux Ratio 
(C>FR). The necessary riser height to ensure the minimum 
recirculation flow required for the CHFR stipulation has been 
coiputed using one dimensional codes.lt is concluded that Natural 
Circulation BWRs of hi^ier powers meets the CHFR criterion and its 
safe c^ceration is definitely possible. 
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CHAPTER 1 


1 


IMTROCXJCnON 

The 1990s* has witnessed the spawning of new nuclear reactor 
designs In US and Japan whlc^ span the technological range from 
evolutionary light water and heavy water, throu^ ll^t t«ater 
passive and gas cooled, to liquid metal fast breeder. 

The likelihood of successful development and operation of aiy 
Individual reactor design depotds on a nix of both ' hard ' end 
'soft ' variables. The hard variables are the more objective 
factors such as actual desiffi characteristics, cost, length of 
technology leap, and degree of research and development support by 
the government and Industrial orgwlsatlons. 

The soft variables are more ambiguous and depend more on so 
called * pid>11c percepticxi * rather than on facts. TTtey may 
Involve Imponderables like attitude towards carbon dioxide 
production and global warming, and alternative power technologies. 
Other factors are more nation specific. Including degree of 
political opposition to siting of nucl^ plwts or radiation 
waste storage sites, pihllc amblvalenpe towards the isricnowns of 
the new technology versus the familiar one and the regulatohr 
procedures. Yet other soft variables are utility specific such as 
size of available marked and suitability of a design's power 
^output* ' 

Mcmt of the designs developed so far have aii^ at Increasiad 
s^ety and ecxhomlcsal; slwllfylri® .cppsti^ 

opera^lbn.-y.. and . 'ma^htBiahce?. becaMse , these' ;dBs1ijy» -"hava. IhclMded. 

/ai$; 



2 


Thf Advanced Bf)131rig W. ttf>r Reactor (ABWfO and its smaller 
c<Hjr»terp.'ir I the? Simplified lioiling Water Reactor (SiBWR) are among 
U)e fvewf^r designs in E?WR proposed by Gerteral Electric (GE ) 
(Reference Nf^ Reactor 16]). l)>e ABWR and SBWR designs are expected 
to strow improvenent over tire current fleet of BWRs in plarrt 
availability, operating capacity factor, safety and reliability 
while reducing power gerreration costs, construction times, 
occupational radiation exposure and radioactive waste. By 
simplifying desi^i of the components system and structure and by 
using natural circulation of coolant further improvements in 
safety, performance and economy can be made. 

The designs of ABWR (fig 1.1) and SBWi (fig 1.2) have much in 
common, the only difference between them are the power rating, 
core flow recirculation (10 internal punps for ABWR and natural 
circulation) and extent to which some of safety systems uses 

I ' 

active versus passive features. An irrportant diaracteristic of 
both ABWR and SBWR is the elimination of external circulation 
piping afid therefore permits a compact containment design. It also 
allows elimination of large vessel nozzles below the core, spd 
therefore the design of a more economic emergency core cooling 
systems. Elimination of external circulation piping results in a 
greater then 50 X reduction in welds, less in service inspection 
(ISI) of the primary system boundary, less occupational dose 
during ISI, increased system integrity and sizable saving in 
capital cost. 

In response to the increasing interest in potential future 
nuclear units cofnprising the characteristics of smaller size, 
greater simplicity and more passive safety features, GE started 
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Advanced 
Boiling Water 
Reactor Assembly 

1 Vcni and Head Sprav 

2 Steam Dr)cr 

3 Steam Outlet Flow Restrictor 

4 Steam Separators 

5 RI^V Siahili/cr 

6 Feedwater Sparger 

7 Shutdown Cooling Outlet 

8 Low' Pressure Flooder (LPFL) 
and Shutdown Cooling Sparger 

9 High Pressure Ckire Flooder 
(HPCF) Spaiger 

10 HPCF Coupling 

1 1 Top Guide 

12 Fuel Assemblies 

13 Core Shroud 

14 Control Rod 
13 Coie Plate 

If) In-C'ore Instrument (iuide 
7 uhes 

17 Ointrol Rod (fuide Tubes 

18 Core DilTereniial Pressure Line 

19 Reactor Internal Pumps (RIP) 

20 Thermal Insulation 

21 C'-onirol Rod Drive Housings 

22 Fine Motion Contiol Rod 
Dnves 

23 RIP Moioi ('asing 

24 LoC4d !V>\\ei Range* Monitor 


Fig. Advanced Bolling Water Reactor 


Assenibly 
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Simplified 
Boiling Water 
Reactor Assembly 

1 Rcacior Pressure Vessel 

2 RPV Top Head 

5 Integral Dr)'cr-Separator Assembly 

4 Main Steam Line Nozzle 

5 Depressurization Valve Nozzle 

6 Chin>ney 

7 Feedwater Inlet Nozzle 

8 Reactor Water CIcanup/Shuidown 
Cooling Suction Nozzle 

9 Isolation Condenser Return Nozzle 

10 Gravity-Drivetr (aM)Iing Ss-stem Inlet 
Nozzle 

11 RPV Support Skirt 

12 Ca>re Top (uiide Plate 
15 Fuel Ass<*njblies 

14 0)1 e Plate 

15 (Control Rod Guide Tubes 
It) Fine Motion Gonirol Rod Drives- 

Slupllfied BorUng Water fteactor 
Assembly. 

f ,y if ' a » 


s 

stilriu''- In 198? of ft 600 BWR with sinir)l1fied powf*r gerjeration , 
safety arid he.a1 rf^irtv/jl systerns. Tt»e basic c^Jtctives that were 
establis1»ed for this new desigr^ SBWR are 
k Power generatiofi cost to be cheaper tlrarr those of coal. 

► Plas'it safety design simpler tlran in current desigri by using 

► passive safety concepts. 

► Design based on existing techrology. 

► Shorter construction schedules. 

► Electrical rating in the 600 MWe range. 

► Inprove the seismic resistance of the core. 

► Lengthen the continuous operating period. 

► Simplify the system by eliminating movable components. 

► Ifnpr~ove the operability and the maintainability. 

Selection of the natural circulation as the means for 
providing coolant flow through the reactor, coupled with a ^2 Kw/1 
core power density, results in a nunber of benefits to help 
satisfying SBWR objectives, Coipared to the existing, forced 
circulaticp plants, the natural circulation BWR offers low fuel 
cycle costs, fewer operational transients, and increased thermal 
margin for the transients expected to occur. In addition, 
elimination of the recirculation, pumps and controls needed for 
forced circulation substantially simplified the design, (fig 1.3) 
Conventional BWR safety / relief valves - wrfiich opens 
closes to discharge reactor vessel steam to suppression pool - are 
not needed in SBWR because in it an isolation condenser is placed 
in the isolation condenser pool so when the reactor vessel is 
isolated frcm the turbine condenser, the isolation condenser 
controls the reactor pressure automatically without the need to 
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of sa*fR. 
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r«ir)ve flu^d f rooi ttif* rfs^r.tor vessel. 

Gravity-dr Ivet) cx)re cooling systerii provide a simrrle approach 
to en>ergBr»c;y core cooling system because it elimirrates the need of 
pumps. Simplif iciation of tf»e cofnpor>ents and the system ma0e the 
primary contaiiyrrent vessel and reactor building snialler which 
reduces the radiation exposure and constructiorr cost 
substantially. 
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CHAPTER 7 

PRINCIPLE OF NUCLEAR REACTCM?S 


2.1 Introduction: 

The thermal energy produced 1n a fission power plant Is the 
kinetic energy of the fission fragments and to a lesser extent of 
the emitted neutron and other particles and radiation SLKsh as 
gamma rays which get converted to heat when these particles are 
absorbed. This heat 1s removed by a coolant and subsequently 
utilized In a thermodynamic cycle. 

Nuclear reactors are variously classified according to 
general purpose or function, type of the moderator, type of the 
coolant, neutron energy classification, type of fuel, type of 
core Internal design and other factors. 

Wille the other possible coolant besides ll^t water cai be 
used like the heavy water, organic liquid, gases and llcyild 
metals, 11 ^t water has be«^ used most extensively , because of 
Its availability, low pumping power and the advanced state of 
knowledge concerning Its chemistry and thermodynamics properties 
etc. 

2.2 Lls^it Water Reactors: 

The reactor 1n which the coolant Is light water may be 
broadly divided Into PWR and BWR. As the names suggest, Bolling 
Water Reactors are those In which the ii^ter bolls within the 
reactor, and In Pressurised Water Reactor the pressure Is high to 
avoid boiling. Thus the water pressure corresponds to the 
saturation tenperature at the reactor pressure In a BWR. The 
reactor pressure Is roughly between 600 and 1000 ps1a. A schematic 
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arrangemerit of a BWR is shown in t^»e figiireC 2.1 ) 

In a Pressurised Water Reaotor , reactor pressures are higher 
and of tf>e order of 2000 psia, for the primary coolant loop within 
the reactor. A heat exchanger in which the primary coolant 
exctianges heat with a secxxidary coolant is used. A schematic 
sketch of a PWR is shown in the figure( 2.2 ) 

However, in this study, we focus our attention on the former 
type of reactor that is the BWR. The boiling reactor has a 
function closely resentling that of boiler in a conventional 
fossil-fuel steam power plant and is basically sin¥)ler than it. 
While in a boiler heat is transmitted from the furnace to the 
water indirectly - partly by radiation, partly by convection and 
partly by conduction: with combustion gases acting as an 
intermediate agent or coolant, in the Boiling Water Reactor, the 
coolant is In direct contact with the heat producing nuclear fuel 
end boils in the same conpartment in which the fuel is located. 

The simplest form of a boiling reactor power plsnt as shown 
in the figure( 2.3 )consists of a reactor, a turbine gaierator, a 
condenser and the associated equipments ( such as ejector, cooling 
system etc.) and a feed pm^ to force the incoming coolant to 
core. Slightly sub-cooled liquid enters the reactor core in the 
bottom wrfiere it received the s^sible heat for saturation plus 
latent heat for vaporisatiai. Vifr®n it reaches the top of the core 
it has been converted to a very wet mixture of liquid and vapor. 
The vapor separates from the liquid, flows to the turbine, dcses 
work and is condensed by the coidenser, and the condensate is then 
punped back to the reactor by the feed punp. That is how the power 
is generated in the conventional BWR. Such a cycle is called 
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Vapor 



'^'<5 i Schematic arrangement 
of boiling-reactor power plant. 
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Hot primary 

coolant Steam 



Feed -water 
heaters 


^ Schematic arrangement of liquid-cooled-reactor power plant. 
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direct cycle BWR, whic^i has the disadvantage of primary coolant 
coming in contact with turbine etc. To overcome ttiis problem one 
can employ indirect cycle BWR in \»^icli a heat exchanger is used to 
separate the rcadioactive primary coolant from the non -radioactive 
secondary coolant. 

2.3 Itetural Circulation BWR: 

The sinplified and modified version of the conventional BWR 
is the natural circulation BWR, more generally knoM^ as Sinplified 
Boiling Water Reactor (^WR). The present thesis deals with this 
sinplified BWl. Except for the neutronic / thermal-hydraulic 
design, there is no difference between the design of natural 
circulation BWR £«d forced circulation BWR. In the natural 
circulation BWR the necessary recirculation flow fore core cooling 
will be produced by the driving force originating from the 
buoysfficy due to density head difference in the closed loop of the 
recirculation flow, (in the Essence of recirculation pumos) 

Actually the technology for the natural circulation is not 
new to BM=5. The Dod-^vaard plant in Netherlands has operated cn 
l?>is principle at a life time capacity factor of 84*1, Larger BWRs 
have been cHperated at 50% power levels in natural circulation mode 
to prove that BW? of this type is indeed possible. 

A schematic of the natural circulation BWR as it contrasts 
with the forced circulation BWR 1s shown in the figureC 2.4 ). 
The major desi^ aspects of the natural circulation BW? has been 
shown in the figure ( 2.5 ). 
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Fig,2.4 Schematic of (a) forced-circulation BWR and (b) 
natural-circulation BWR. 


Coro Miser 
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Flgire (2.5): Major design aspects of natural circulation BWR 
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CHAPTER 3 

PROBLEM STATEICNT A^D PROCEDURE 

For deterniining the feasible power rates for the 
(natural circulation BWR , the design criteria for steady thermal 
hydraulics 1e the critical power and core stability should be 
satisfied. The feasible power rates from the point of view of 
critical power requirements has to be determined. 

FLOW MOCEL: The flow In the riser .downcomer and 
parallel core channels was assuned to be vertically one 
dimensional. A code was developed to find out the minimum core 
flow necessary to meet the critical heat flux requirements. The 
Input parameters of the code are the total power. Inlet flow sub 
cooling and the core operating pressure. 

Next a code was developed i»(h1ch cormxites the total 
pressure drcH^s In the respective parts of the recirculation flow 
formed by the dowHxximer -core & riser and the necessary riser 
he1^t to ensure the driving pressure for natural circulation. 

The critical heat flux ratio has been taken to be 


greater than 1.3 from safety considerations. It 1s assumed that 
the selected requirements for CHFR are conservative «xxigh to give 
an over all Inherently safe reactor . 



tl»p proos'Ckjre hfj^. been outlined in the following fig 



figure S.lsoutHne of procedure for 
calculation of riser height 
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CHAPTER A 

CRITICAL HEAT FLUX REOUIREMEhlT 

4.1 Introduction 

Tfie liquid coolant used In nuclear reactors Is 
sub 3 ected to h1^ rates of heat fluxes .In boiling water reactors , 
the change of prfiase Is deliberately designed Into the system atf^d 
should be controlled. The ccx^dltlcxis at which the heat transfer 
coefficient of the two phase flow substantially deteriorates 1s 
called tine critical heat flux condition (CHF) For a given flow 
conditions, 1t occurs at sufficiently h1^ heat flux or wall 
temperature. In a system, in i»^1ch the heat flux Is Independently 
controlled , the consequences of OF occurrence 1s the rapid rise 
1n the wall temperature .For systems In which the wall tenperature 

I, 

1s controlled , the occurrence of OF Inplles a rapid Jifxjrease In 
the heat flux. 

Many terms have been used to denote the OF 
conditions. Including "boiling crisis" and " burnout". Departure 
from nucleate boiling (DJ®) Is the term originally used to 
describe the OF conditions In pool boiling. It can also be 
encountered In flow boiling , when the formatlcxi Is rapid enough 
to cause continuous vapour film at the wall. Generally the OF has 
been calculated using both the criteria of (a) Departure from 
nucleate boiling aid also from (b)dryout. 

4.2 Critical Power Ratio(CPR) 

To eliminate the undesirable features InheoKit In 
the local CHF hypothesis a new correlaticwi known as the General 
Electric critical quality- boiling length (GEXL) correlation was 
developed.The GEXL correlation Is based on a large amount of 
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bDillng transition data taicen in General Electric's ATLAS Heat 
Transfer Facility .The Generic form of the GEXL correlation Is 

Xc.r = Xct( lb, dh,Q,l,p,R) (4,1) 

where xcr = bundle averaged critical quality 
lb = boiling length 
dh = equivalent diameter of channel 
1 = total heated length 

p = system pressure . 

I, ' 

R= parameter that characterises the local peaking pattern 

The GEXL correlation Is a "best fit " to the 
experimental data and said to be afcle to predict a large range of 
data. The GEXL correlation Is a relation betwe«> parameters that 
depend on the total heat Infxit from the channel Inlet to a 
position within the channel. This correlation line Is plotted in 
terms of critical bundle radially average quality vs the boiling 
length- For a g1v«i bundle power ( Oi.or Q2) the locus of bundle 
conditions can be represented in the fig 4.1 The critical power of 
a bundle Is the value of bundle power to which the correlation 
line becomes becomes a tangent. The critical quality - boiling 


length leads to a Identif 1catic« of a new margin of safety In 


designs 1e the Critical Power Ratio. CPR = 


critical power 
operating power* 


4.3 Critical Heat Flux Ratio Correlation 


In abs«x» of the GEXb correlation ,(1t being a proprietary 


Information) we take recourse to the OF Correlations as given by 


Biasi.dodreas & Kazimi t?]) 

-l/Cr “1/6 -I 

qbIa8i=(2.?64-10^)(100D)'"G h.468F(Pj^^^)G - x Jw/m 

....(4.2) 

qbia»i=(15.O48-10’^)(lO0D)‘V^'®H(pbar)n- xl W^m^ (4.3) 


I 



tangent point 



Figure ^ Correlation of criii 
conditions under dryout. 


The greater of the above two values is taJcen where 

F(pbar) = 0.7249 +0 .099 ' (pba r ) -expC” 0.032 pba r ) . . , (4. 4) 
H(pbar )=-l . 159 + 0.149 Pbar(exp(-0.019 pbar) 

+ 9 pbardO + pbir)'^ ...(4.5) 
pbar = lOp where p Is In Mpa 

the equivalent diameter D = 4*area/ wetted perimeter. 


n=-{0.4 wrfien D^fl.OI 
0.6 when EK0.01 

The quality x(2) has been cxanputed as 

qcr A heub 


x(z) 


...(4.6) 


mhfg hfg 
where x(2) = quality at any hei^t 2 

til = mass flow rate 

Ah sub * degree of sub cooling 

The pattern of heating has been taicen to be uniform and the axial 
peaking factor to account for the sinusoidal heat addition has 
been taJ<en to be 1.5?.An iterative scheme was developed as under. 


assume a total core flow G 


:e 


G = G + AG 


c::j#J.CULATION OF CHFR by Biasi relation 




min CFFR >1.3 


IF NO 


IF YES proceed 

fig 4.2:outline for calculation of min G 


For the two power densities of 42 and 50 kw/lit the (iHFR is 
plotted against G . Also the minimum 6 retjiired for CHFR to be 
greater than 1.3 for different fuel lengths is shown in the 


Table 7.2 




CHAPTER 5 

RISER HEIGHT CALCULATION 
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Wtnle Uie last c^^apter looked at the necessary 
minTinum core flow to ensure that the reactor was safe from the 
critical heat flux considerations, the present chapter evaluates 
the r>ecessary riser height to ensure the required core flow by 
natural circulation .The equations developed in this chapter are 
from El-Wakil. [53. To conrpute the necessary hei(^)t the approach 
has been to equate the two s 

(1) pressure drop components as the coolant flows through the 
various components and 

(2) the gravity hea3 that is available because of the density 
differential between the liquid in the core& riser and the licyiid 
in the downcomer. 

5.1. ratio of non - boiling to boiling hei(#it 

In order to correctly evaluate the pressure drops 
as well as the average density in a boiling channel ( necessary 
for evaluating the driving head in a natural circulation reactor ) 
it is required to calculate the boiling and the non - boiling 
heists in the channel. The non - boiling hei^t Ho is that in 
which only the swsible heat is adcted to the incoming sub cooled 
ccolant at the channel bottom. The remainder of the channel is that 
In which the boiling tskes. place. ( hb). 

The ratio Ho/H (H: total Iwigth of channel) may be 
evaluated from the ratio of sensible added to the total heat 

added.as given below. — — KL . 7 .. lli — (5.1) 

” (hf + X • hfg) - hi ^ 

e 

where qs = sensible heat added per mass flow rate of incoming 
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coc)lar»t 

qt = tx)tal htvit adckni per mass flow rate of incxamlng fluid, 
hi = enthalpy of 1r»can1ng coolaf>t to thie core. 

In case of sinusoidal heat addition (extrapolation lengths 
neglected) the area bounded by z = Oboz = Ho 1s proportional 
to the Qs and ttie total height to the qi.It can be written as 

Ji Ho. 


qg - 
qt 


q^ s1n( ) dz 


j” q^ s1n( ) dz 


= 0.5-( 1 - cos( 


H 




(5.2) 


5.2 pressure drop components: 

5.2.1 friction drop In core channel 

In the boiling hel^t Ho , friction Is due to single phasedlquid 


only) flow and may be evaluated by Darcy's formula. 

Ap«P=f ^ (5.3) 

where /)psp = single phase pressure drop 

f = frlctloTi factor dependent upon Reynold's no. 

De = equivalent dia of diannel given by 
4" diannel area/(wetted perimeter) 
p = average density In non - boiling height 
v= average velocity 1n the non - boiling height 
pressure cfrcp in boHing height : 

In the boiling hel^t the pressure drop due to friction Is greater 
tfw) that die to single - phase liquid flow for the same flow 
rate. 

The two phase drop as It compares with a calculated 
friction drop, assuming that only saturated llc^ild existed In the 
chaainel Is given by R = — - •• -(5.4) 

where R = two phase friction multiplier 

0 

iiptp = two phase pressure drop In boiling height only 
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= saturatf»d l^cpild phase pressure drop In the saine 
tteigtit 

R Is found to be a function only of the ©xH void 
fiar^tion and independeiit of the many operating conditions such as 
pressure , ^leat flux , coolant flow rate and slip ratio. 

The Lottes' Flinn correlation has been used which is given by 



wt>ere a® = exit void fraction 
5.2.2 acceleration pressure drop 


As the coolant receives heat in the core channel , 
it undergoes an increase in volume, ccxisequently it has to 
accelerate as it travels through the channel .The acceleration 
pressure drop has to be evaluated from Ihe relationship involving 
the change in momentum of the incoming a^d outgoing fluids. Force F 
due to the change in the momentum Is given by 
F = Apa Ac= Ifif Vf® + fhg Vge - iflt VI (5.6) 

where Apa = pressure drop ckje to acceleration or change in 
momentum, 

Ac = crossectional channel area, 
vi= inlet velocity to channel 

Vfe, vge, = velocities of saturated liquid and of vapour 


respectively 

enyj 5.6. can be derived as Apa = r*G^where r is 



0 


= saturated specific volune 
vg = saturated vapair velocity 
vf = saturated lic^id velocity 
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5.2.3 expansion losses in two pftase flow 

Eie4>=itise of sucJdw* area ctianges as they flow in 
reactor ct«ir>els, two rfjase mixtures undergo cftanges in pressure, 
for eg. w^>eri stK.-h mixture leaves the fuel channel to enter the 
riser . (comti«ori to all tfre fuel elejiients) . Force balance across the 
expansion is written as. : (fig 5.1) 

pi A2 + mfl-Vfl + mgl-Vgl= p2-A2 + mf2’Vf2 + mg2-Vg2 

the eqn reduces to 

r 1 1)2 

= i - a^ay -~rr 

A5 

(5.8) 


a) 


P9 




5.2.4 contraction losses in single phase flow 

The area ccxitraction results in a vena contracta of 
area Ao following the oontracticxv (fig 5.2). C^ontraction loss is 
assunred solely due to the expansion of the fluid from the vena 
contracta to tire point where the fluid fills the area A2.The 
contractiw loss occurs when the fluid enters the fuel 
chanr^ls.Cfrom the bottom) For single fhase it may be written as 

P2 - pi = 0.7-f--^ 

^ A? A^^ 

where iht is the mass flow rate . i > . iti 

5.3 DRIVING PRESSURE IN NAIURAL CIRCULATION - ' « ' 

The water in the downcomer is either saturated or sli^tly sub 
cooled, cMaending upon where the feed water is added to 
recirculation water -Let the density in the downcomer be given by 
p The riser hei^^t is ai unheated extension of the core but 
usually with fewer walls, dividers, and ofcourse , with no fuel 
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Fig. . 


Pi 



Single phase flow In sudden contraction 
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s so that t.tw*tf* H 3f‘ss fricticxi thatfi In the core .The 
qij.lHy In ttw riser is substantially Ure sarne as that at the 
ctururt+l exit, tint is x = x Tt»e slip ratio in the riser is 
assufiied to be sanre as in the core .Ttie density along the riser is 
als(j tire saifie as ttrat at the core outlet. It can be readily seen 
tJrat ttre addition of riser also increases the driving pr^essure by 
the quantity g.The driving pressure in the case of 

ctrannel with chirmey is given by 

^ ^dc^” + pT] g...(5.10) 

where p= average density in chanr^l given by 

(p^Ho + Hb)/ H where p^ is the average density In the 

non - boiling charrnel,p^ the average density in the boiling 

cfiainel. p^ may be closely evaluated as the inverse of the 

average specific volume in the non - boiling channel as follows 
— =0-5* (v ,+v,). 

P I ( 

O 

For the case of sinusoidal heat addition in the 
channel an expression of p^can be derived .Taking a differential 
elemerrt of height dz at hel^t z within the boiling region ( that 
is 2 > Ho) The defrsity of the two phase mixture is given by 
( 1- a^) P^+ «^Pg = Pf- a^( Pf Pg) ..-(5-11) 
average density will therefore be given by 

To integrate this equation , a relationship between (the 
void fraction at z) and z roust be obtained - 

.First a relaticxiship between and x^.the 

X 

PMality at z is obtained as ^ --(5.13) 
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V-, ricMtf' by by 7 tn hf-il addition as follows 

q.' = Chf 4 x^-hf<i> - hi ..,{5.1A) 
qi - (hr 4 X flic,)- hi ..(S. 15) 
wti^'re q/ = Not, gt'wr/itpd and transferred up to hel^t z and 


qt = total generated (per mass flow rate). For sinusoidal heat 

addition , the Oi'>e under consideration: the above two quantities 
are related by eet expression as — ^=0.5-( 1 - cos — ^ . (5. 16) 


qt 

The solution to Integral Is obta1r>ed as follows 
(a) for 


^ ' Or 


_i! . ^i'^4 ~ ^ 2 ‘S H 

T ' ^ 


' 4 J i i~ 
’ c - c 


HB 


H 




^ tan-’ 


n 




(b) for > Cg 

= p- ip. -py 

hi * 9 


2 2 

C - C 
3 4 


C *0 - C *0 

1 4 2 3 H 


2 

''4 


HB 


In 


( c - c,) tan( 

4 3 


«H0 

2H 


)* 1 - 


where c 


^ 2 = - 


3L 


( c ~ c,) tan( 

4 J 

hf - hs 


TtHO 

2H 


) - 


2 

^ - 


1 2hfs 
qt 


hfs 


2 hfg 

Cg= ip 4 (1-tp) 

c = (1- tP)c, 


....(5.19) 


j...(5. 


18) 


The density of the liquid in the downcomer is assumed to be the 
same as that of the saturated liquid becaise of the h1^ 


recirculation ratio- 
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f.HAf'TtR 6 

FLUID Flow PATTERN IN MULTI CHANNEL CORE 

III ctw 4 »t.t*r S the riser heiglit has been cjon^ted 
asMJitirio that lt»e hncit gerjerated In the core is uniform 
radially. Ac;tu<'illy , the neutron flux distribution in the radial 
directicxi is r>ever re.illy uniform; the amount of heat generated by 
the the fuel differs from orie chanriel to another. The amount of 
stec«» generated and consequently the exit c^lity vary from one 
channel to arrother. In a natural circulation reactor, the driving 
pressure is therefore not equal in all channels, being greatest 
for those cLvannels with the highest neutron flux( ie near the core 
center) arid least near the core periphery, The coolaait flows in the 
different channels also vary in the same order. 

For calculation of the fluid flow pattern , the 
core is divided into a no. of regions, each containing a no. of 
fuel channels. These region are chosen so that each contains fuel 
channel at the same neutron flux level.These regions do not 
necessarily contain the same no. of fuel channels, nor are they 
necessarily adjacent to each other, for exauple in the fig 6.1 a 
72 fuel assembly has been divided into 4 regions containing the 
following no. of fuel channel 

region no of ft^l channel 
A 16 

B 20 

C 24 

D 12 

In our case the core hats been divided into 7 regions . A code was 
developed that computed the distance of the centroid of the fuel 




char.wl front th« core c;ertter. and accordingly divided the core as 


RANGEdn m) NO (T FUFL CHANNELS 


A 

0 - 

0.6 

60 

B 

0.6 

- 1.0 

104 

C 

1.3 

- 1.6 

112 

D 

1.6 

- 1.9 

176 

E 

1.9 

- 2.1 

224 

F 

2. 1 

- 2.3 

100 

G 

2.3 

and above 

96 



Total = 

872 


For calculation 1t Is assumed that the data given are the-fehe- 
total core thermal power , the cc^ (^jerating pressure, the inlet 
flow sub cooling, the general core design .aamd the core shape.The 


quantities to be computed are the riser height i^d the core flow 
pattern. The steps outlined are 

(1) The heat generated in each of the regions has been calculated 

using the total core power o«d the neutron flux distribution. The 
heat ger»erated is assigned to vary according to Jo (the Bessel 
furction). i, '' ' 'f 

(2) the core average exit quality has been computed from the total 
core thermal power and assuming the total flow through the core. 

(3) Starting with the region of max. heat generation , a value of 
the exit c^iality is assumed . 

(4) The total flow In the hot regies is calculated using the total 
heat gerierated . Using the axial flux distribution (sinusoidal) 
the non -• boiling and tte boiling heights are convxited.The 
pressure losses in that region are computed. 

(5) TL« pressure drops in the downcomer are computed using the 
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lot<»l flow. 

(6) T^li• riser Is cxxnr^Jted using Ue procedure outllt'ied in 

clrafiter S 

(7) An exit cjiiality for otter regions is assumed. The flow af)d 
cx>r respcxrding pressure losses in that region is conxauted using the 
riser hei^t in (6). From the xe assimed arid the axial flux 
distribution , the average density , and the driving pressure for 
that region is calculated. This driving pressure is conpared with 
the losses in that region plus those of the downcomer. If the two 
values are not equal then it Is necessary to assume eawther value 
for X® ar»d recalculate until equilibrium is attained. 

(8) This step is repe^ated for eacdi of the regions. This fixes the 
value of flow ^d tte exit quality at each of the regions. 

(9) With these values the total core flow and the average exit 
quality are calculated. This is then compared with (2) 

Steps (2) to (9) are repeated till the two values are ec^ial.This 
gives the core flow pattern arid the riser hei^t.This procedure 
has been done for two values of core average flow of 6.0 & 9.6 
kg/s bundle. The flow distribution has been sho«i in Tables 7,5 & 
7.6. 



CHW>TER 7 

RtaJl TS AND DISCUSSION 

7.1.1 varlatiioo of CHFR with core flow 

As Uie ref;irculat1on flow Increases , the critical 
flux ratio is fcxifid to increase .This is sho^n in fig 
7.1,w!»erein tlie CHFR increases with the recirculation flow. This is 
dcxx' for two pcwer densities of 42 & 50 kw/lit. Therefore the 
n»atgin of safety according to CHFR Improves with increase in 
recirculation flow. The chaaage in minimun recirculation flow with 
change in active fuel lengtlj has been shown in the fig 7.2 (for 
the CHFR to be greater thatf^ i.3) 

7.1.2 necessary riser heig^it 

Following the proceckire outlined in the chapter 5, 

tlie riser heigfit was computed to produce the necessary minimum 

recirculation flow This has been done for the power densities of 

42 & M kw/lit These are shown in the fig 7.3 .The riser height is 

found to increase with the active fuel l^^gth but at hi^r power 

densities it irtcreases at a much more rapid rate. The change in 

riser hei^it with power densities is shown in the fig 7.4.This is 

done for different active fuel lengths . According to the results 

of this study, the necessary riser hei^t is mainly dependent on 

the active ftj^l length(L) aid the core power (cj) 

.The riser hei^t «f(L,q) 

An increase in power rates is possible with increase in riser 
height. Yet the riser hei^t cannot be Indefinitely increased 
because of economical considerations -Thus the riser height has 


been limited to 10m 
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TAfy F 7.1 :CHFR VS recirculetlot) flow 


c, 

(Kg/s._biif>rilp) 

pd ' 4?kw/lU 

pd * 50kw/lit 

(n R 

(PR 

CPFR 

CPR 


result 

result 

In refl 

result 

obtained 

result 
in refl 

A. 000 

0.739 

1.0 

0.621 

0.85 

5.000 

0.891 

1,12 

0.748 

0.95 

6.000 

1.0?7 

1.24 

0.863 

1.05 

7.000 

1.148 

1.35 

0.964 

1.15 

8.000 

K?5S 

1.45 

1.054 

1.20 

9.000 

1.349 

— 

1.133 

1.25 

10-00 

— 


1-201 

1.35 

11.0 

— 

— 

1.261 

— 

1?.0 

— 

— 

1.311 
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tiJhO 


6(kg^s*bundJeJ> 

CHFR VS MIN RECIRCULRTION FLOW 
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TAPIF 7.7 : ni’idlfmiiii G vs active fuel length 


ar.t We fuel 
leiKjthd) 

minimum flow rate 

pd- 42 kw/lH 

pd = 50kw/lit 

( In ni) 

kg/s bundle 

kg/s bundle 

7.40 

7. SO 

10.10 

7.50 

7.90 

10.70 

2.60 

8.20 

11.30 

2.70 

8.60 

11.90 

2.80 

9.00 

12.60 

2.90 

9.40 

13.30 

3,00 

9.80 

14.00 

3. 10 

10.20 

14.80 

3.20 

10.60 

IS, 60 
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TABLE 7.3 : riser hel^t vs aotive fuel length 


active fuel 
lengthCln m) 


2.40 

2.50 

2.60 

2.70 

2.80 

2.90 

3.00 

3.10 

3.20 


RISER HEIGHT (In m) 


pd=42kw/11t 

result 

obtained 

from 
ref 1 

result 

obtained 

from 
ref 1 

3.59 

3.80 

1.38 

1.33 

4.32 

4.50 

1.76 

1.67 

5.12 

5.50 

2.li 

2.0 

5.99 

6.40 

2.58 

2.5 

7.04 

7.50 

3.08 

3.0 

8.19 

9.0 

3.61 

3.4 

9.45 

10.2 

j 

4.19 

4.0 i 

10.94 

11.20 

4.82 

4.6 

12.60 

— 

5.85 

— 









RISER HEIGHT VS RCTIVE FUEL LENGTH 
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TABLE 7.4: riser hel^t at different power densities 


power 

density 

(kw/llt) 

RISER LCIGHTdn m) i 

! 1= 2.4m 

1 

1=2. 7m 

1 = 3-201 

42 

1.38 

2.58 

5.85 

i 

43 

1.61 

2.88 

6.1 

44 

1.81 

3.25 

6.74 

i 

45 

2.05 

3.64 

7.41 

46 

2.35 

4.04 

1 

8.22 

47 

2.62 

4.46 

9.07 

48 

1 

2.89 

4.9 

10.08 


3.23 

5.43 1 

11.15 

50 ! 

3.59 

5.99 

12.6 
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7.1.3 flow distribution In the core 

The flow dlstr1butl<x> In the core channel has been 
calculated according to the procedure outlined In Chapter 6. This 
has been dorte for two total core flows of 6 end 9.6 kg/s 
bundle. The flow distribution obtained has been shown In the tables 
7.5 & 7. 6. The riser height confuted has also been shown In the 
tables.lt clearly shows that the core flow will be the highest In 
the core channels with the highest heat generation (corresponding 
to the core channel with the hicfiest neutron flux distribution.) 

7.2 CONaUSION 

The thermal hydraulic requirement for CHFR was 
studied for the natural circulation boiling water reactor. The 
necessary riser height was calculated to establish the required 
flow In suc^i reactors.lt Is concluded that natural circulation 
BWRs of higher power rates will be feasible even If major 
restrictions on the reactor pressure vessel diameter and riser 
hel^t are set .Also the results cire In general agreement with the 
published results of Akira Yasuo et al.(l) (refer to table 7.2 for 
comparison with published results) 

Tte main caicluslons are as follows 

(1) With increase in core flow, the critical heat flux ratio 
Increases. 

(2) The core flow carj be maintained using a certain riser 
height. 

(3) The necessary riser height 1s a simple function of power 
density asid active fuel length. 



TABLE 7.5 :flow distribution pattern 


for average flow of 6.0 kg/s. bundle 


REGION 

FLOWdn kg/s. bundle) 

A 

8-74 

B 

7.83 

C 

6,90 

D 

6.14 

E 

5,39 

F 

4,36 

6 

3,32 


valu of riser he1^t 1s « 2.5in 
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TARl E 7.6:flow distr ^but1or> pattp'rn 
for average flow of 9.6 kg/s. bundle 


REGION 

FLOW 

A 

15.79 

B 

12.68 

C 

10.85 

D 

9.5? 

E 

8.27 

F 

6.63 

G 

5.02 


valu of riser height 1s 6.5m 



CHAPTER 8 


RECOWCNQATIONS FOR FUTURE WORK 


(1) It Is desirable to use the GEXL correlation for determining 
ttie Critical Power Ratio (CPR) as 1t will give more accurate and 
dependable results. (CPR Is defined In chapter 4) 

(2) Tte model taken for analysis In this study Is one - 
dimensional. The flew distribution should be calculated using a the 
radial heat distribution (for eg. the Bessel function) while 
satisfying the Critical Power Ratio requirement. The riser hei^t 
should be calculated accordingly as outlined In Chapter 6. 

(3) The heat distribution 1n the axial direction has been ta^en to 
be strictly sInusoldal.The work can be extaided to other heat 
distributions which are not strictly sInusoldaKas Is the case 
actually) like the c^Kspped cosine function, and a heat distribution 
in which the maximum value has been shifted closer to the bottom 
as in actual reactors. 

(4) The pressure drops should be calculated more extensively by 
taicing into account the spacers, support plate loss coefficient; 
also examining the case of riser height with separators. The 
accuracy of the riser hel^t calculatiw is affected greatly by 
the accuracy of the pressure drop calculation. 

(5) The study of change of parameters like pressure, no. of bundles 
should be done to get a more detailed idea of such reactors. 

(6) The case of Natural Circulation BWRs with external 
recirculation can also be examined. 


operating pmssur©(MPa) 

7.0 

core Inlet sub cool1ng(KJ/kg) 

33.0 

core diameter 

5.8m 

fuel bundle type 

8x8 

no of fuel bundles 

872 

active fuel length 

2.4-3.2in 

radial peaking factor 

0.65-1.4 

axial peaking factor 

1.57. 

fuel pin diameter 

14.2nin 

fuel channel dimension 

0- 138x0. 138ra 
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APPENT^IX B 

{ This prf>gr<»ii J'lar-. bef;*rt used tx» calculate the riser hels^it.The 

ffilnltmui rec:1rculat1w has been conyiuted by the prcx::»cfcjre crIt.The 

Functlori riHi Is used to corffwte the non- bollino height. The 

procedures, btidl and bhd2 confute the average boiling densities In 

the core. The riser height has be^ finally calculated by equating 

the pressure drops and the driving pressure.} 
program natur_c1 rc ( 1 nput .output ) ; 
const 

p1=3. 1415927; 

hf=1258-18; 

hg=2767.16; 

hfg= 151 1.48; 

hfeecN168; 

vf*0.0014; 

v1=0.0013; 

s=3.5; 

pf =7 14.28; 

pg=35.71 ; 

nob=872; 

mul=2.48E-03; 

var 

cl .02 ,c3 ,c4: real ; 
d: real ; 
deltah:real; 
ct)r:str1ng; 
g.pd.lrreal; 

vfg.vg.hi ,m,mg,qs,qt:real; 
f 1 ,f2,ho,hb:real; 

1: Integer; 

^n1t,A1 ,a2,alphae,vl:real; 

reyn ,prf ,f rf ,R_, rac .accprd: real : 

hyd:real; 

sppad .tpprd: real ; 

pe,s1 ,pb,po,xe;real; 

frch,frdc:str1ng; 

remprd ,dco ,pave ,hch : real ; 
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fut«::i\w> nlitKpi ,l,cr5,ql:newil):re,'jl; 

var 

x,y,u,v:rajl : 
begui 
x:=cjs/qt ; 
y:*l-2**x; 
ui=sqrt(Hy*y); 
v:=u/y ; 

nW» : =arctan (v )*l/p1 
end; 

function bhd1(p1,pf ,pg.cUc2,c3,c4,ho,hb,l:real):real; 

• var x,y,z:real; 
begin 

x:=s 1 n((p 1 *‘ho)/( 2 -{ri))/cos((p 1 *ho)/( 2 .tri)>; 
y -. = In ( (c4-c3 )*x +sqrt (c4*o4-c3*c3 ) ) 

- In ( (c4-c3 )"‘x-sqrt (c4’^34-c3’^:;3 ) ) ; 

Z:=y*'((cl’*l24-c2’^^)’*‘l)/(c4*sqrt(c4*'c4-c3*c3)*hlf‘p1 ) ; 
bhdl : =pf “ (pf-pg)’*'(o2/c4-2 ) 


end; 



prxx»dure bhd2(pb,p1 ,pf ,pg,c1 ,c2,c3,c4,ho,hb,l--«^l)s 

var x,y, 2 :real; 

begin 

x: =s1n(p1*ho/C2 .0*1) )/cos(p1’*ho/(2. 0*1) ) ; 
y: = 1 .0-2.0/p1*arctan((c3-c4)*x/sqrt(c3*c3-c4*c4)) ; 
z : = (c1*c4-c2*c3 )*l/(c4*sqrt(c3*c3-c4*ic4)*hb) : 
pb:=pf-(pf-pg)*(c2/c4+y*z) ; 


end; 



procedure cr It (var pd,l,g:r6al); 

const 

pb=70; 

al *0.0088; 

d=0,010; 

hs*33; 

hfg=1511.48; 

var 


hbireal; 



07,qt f .cpr ; 

Uitegf^r ; 
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bf'gifi 
qpr: =0.0; 

Q:=1.0; 

while (cpr<=1.30) do 
begin 
g?:=g/a1; 

hb:=-l. 159+0. 149*pb*'exp(-0.019*pb)+9’*pb/(10+pb*pb) : 

2 : = 1 . 0+ 4* 1* 1 5 . 048*exp ( 7* In ( 1 0 ) I’^xp ( - 0 . 6*ln ( g2 ) )*hb/( gZ^cThf ^ 1000) 
cpr : = 150 . 48’*fexp( -0 . 6* In ( g2 ) )’^hb* ( 1 +hs Aii'g ) /z ? 

Cpr : =qcr/( pd/42000*0 . 802 ) ; 
g-.=g+0.1 
end 
end: 

{ _} 

b6g1n{ main program} 
pd: =42000: 

while Cpd<=50000) do 
begin 

wrltelnCfor power density' ,pd): 

1:=2.4: 

while (1<=3.2) do 
begin 

cr1t(pd,l,g): 
cfco:=5.8; 
vg:=l/pg; 
m:=g*nob: 
deltah:»33-0; 
hi :=hf-deltah: 
d:=(h1-hfeed)/(hf-hfeed) : 
mg: = (pd*p1’*dco*dcxy''l/4. 0-nf*'deltah ) /hfg ; 
xe:=mg/m; 
po:=2/(v1+vf ); 
qt:=hf+(xe*hfg)-h1 ; 
qs:=delta4i: 

cl :=qt/(2*hfg)- (hf-hl )/hfg: 
c2:=-qt/(2'*1ifg); 



hyd:*O.Oin?^i; 

Al:--n. 01907: 

a? 0.008877: 

giiriit:=o/a7; 

vl:=0.9S*g/6.0; 

rey f t : ■= hy cT pf* V 1 ; 

alpitae: =xe/<sH ( 1 .0-s1 )*xe) ; 

pe : = a lp^^ote* pg+ ( 1 . 0-alphae )*pf ; 

f rf : = 0 . 1 8V©xp ( 0 . 2* In ( rey n ) ) : 

R_: = 1 .0/3. 0*( 1 .0+1 .0/( 1 .0-alphae)+sqr ( 1 .0/( 1 .0-alfrfiae) ) ) ; 
rac:=vf*'<sqr( 1 .0-xe)/( 1 .0-alphae)+sqr(xe)/alphae*vg/vf- 1.0); 
accprd:=rac*sqr(gunit) ; 

sf^3cd:=-0.7*( 1 .0/A1*A1-1 .0/(a2*a2) )*g/pf*g: 

tpprd:= ( 1 .0/(Al*a2)- 1 .0/(A1*A1 ) )*g*gi*(sqr( 1 .0-xe)/(pf*'( 1 .0-alphae) )+xe’*‘xey 

c3:=s1+(1.0-si)*c1; 
c4:=(1-s1)=*‘c2; 
ho:= nbh(p1 ,l,qs,qt); 
hb:=l-ho; 

prf:=fr-P*'pf^sqr(vl)*(ho+RJ^hb)/(2.0*hyd): 

If sqr(cA-)>sqr(c3) 
then begin 

pb:=bhd1 (pi ,pf ,pg,c1 ,c2,c3,c4,ho,hb,l) : 
end 

else begin 

bhd2(pb,pi ,pf ,pg,c1 ,c2,c3,c4,ho,hb,l) ; 
end; 

f 1 : =f rf^sqr (in*4 . 0/( pi*p6*dcxy^dco ) )*pe/ (2 . O^dco ) ; 
f 2 : =f rf*sqr (nf'4 . 0/ ( pf*pi=*' 1 0.456) )*pf / (2 . (Tl . 056 ) ; 

renprd:=accprd+sppad+tR3rd+prf : 

pave: = 1 /l’^(pd’'t»o+pb*hb) : 

hch : = (9 . 8 1*pf*l-9 - 8 1*pave*l-f 2*1- remprd)/(f 1 +f 2-9 . 8 1*pf +pe*9 . 81); 
writelJiC M:10:2,' ■,g:10:2,’ hch: 10:2); 

1:=1+0.1 

end; 

pd:=pd+1000 

end: 

end. 
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Af'JlNTJlX C 

{Th^. program dW1(**s Uh=* core 1r»to deferent regions 

^wid rx>rjnt.s rto. of grids In different regions} 

program gr id_crit(fXJtput) ; 
const. 
dLx=0. 138; 
type 

string=array[ 1 . -2193 of real; 
var 

cx,cy,cd:strir>g: 

1 ,3, k: integer; 

cntO ,cnt 1 ,cnt2 ,cnt3 ,cr»t4 ,cnt5 ,cnt6 : i nteger : 

^ } 

begin 

CX[13:=1.863; 
cy[n:=0.069; 
for i:=2 to 196 do 
begin 

3:=(i-1) mod 14; 
cx [ i 3 : =cx [ 1 3 - 3*dtx ; 
k:=i div 15; 

cy[i 3:=cytl 3+k*dtx ^ 

aid; 

for i:=197 to 207 do 
begin 

cx[1973:=cx[13-3’^dtx; 

cyCi3:=cyCl833+dtx; 

3-.=1-197; 

cxt11:=cx[1973-3*dtx 

end; 

for i;=208 to 218 do 
begin 

5-. = 197+i-208; 
cxC13--=cx[33; 
cytil:=cy[33; 
end: 

cnt0:=0; 

cntl:=0; 
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crit ?: -0 ; 
t;nt3:-n; 
t:nlA:^0: 
cnt5:*0; 
c;rit6:=0; 

for “i: = l to ?18 do 
begirt 

cd[ 1 3 : =sqrt(cx[ 1 ]’*cx[ 1 ]+cyC 1 3*cyE 1 3 ) ; 

If (cdt13 <0.6) then 
cnt0:=cnt0+l ; 

If ((0.6<=cd[n)and (cdCil < KOO)) then 
cnt1:=cntl+1 ; 

1f{( 1 . 00<=cd£ 1 ] )and (cd[i3 < 1.30)) then 
cnt2:=cnt2+ 1 ; 

1f({ I,30<=cd[l3)af4d(cd[13 < 1.60)) then 
cnt3:=cnt3+l; 

if(( l.60<=cd[11)and(cdt13 < 1.90)) then 
cnt4:=cnt4+l ; 

1f(( 1.90<=cdn])and(cd[13 < 2.10)) then 
cait5:*cnt5+l ; 

if(( 2.10<=cd[l3)and(cd[1] < 2.40)) then 

cnt6:=cnt6+l ; 

end; 

writelnCthe valu of cntO 1s',cnt0); 
wr1teln(’the valu of cntl Is ',cntl): 
writelnCthe valu of cnt2 is ‘,cnt2); 
writelnCthe valu of cnt3 Is *,cnt3); 
wr1teln(‘the valu of cnt4 is *,cnt4); 
writeln(*tl^ valu of cntS is ' ,cnt5): 
writelnCthe valu of cnt6 Is ’,cnt6); 
end. 
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A('E'm>!X D 

{This H for tifiictil'it kxi of fluid flow distribution 

p. it torn } 

ptogi am rvitur_circ(iriput,out.pit) : 
type 

lisn=ar r ciy[ i. .7] of real; 

var 

k: Integer: 

pof ,nob,xe,inl ,hch,g1 : listl ; 
roav ,xav ,mcure .xclbii: real ; 
procedure rsrht(var xe, pof ,nob,ml ,hch: list!) : 
const 

pi =3- 1415927: 
hf=1258.18; 
hg=2767.16: 
hfg= 151 1.48; 
hfeed=168; 

1 =2.7: 

vf=0.0014; 

v1=0.0013; 

s=3. 1 ; 

pf=?14.28: 

pg=35.71 ; 

qv=42000; 

niul=2.48E-03: 

var 

cl ,c2 ,c3 ,c4: real ; 
d: real ; 
deltah:real: 
cfeirireal; 
pd,g:real; 

vfg,vg,hi ,nv,nig,qs,qt:r6al: 
f 1 ,f2,ho,hb:real; 

1: integer: 

gun1t,Al ,a2 .alj^iael ,alph^,vl:real; 

reyn , prf ,f rf ,R_, rac ,accprd: real : 

hyd:real: 

sppad , tpprd: real : 
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vl 1 ,f f f J ,r«y»t1 ;re/il ; 

ff ch,f rH(;:r«Vil ; 
rtsijftj d: rfvsl : 
totpr d: r«il ; 
ckxt: retd ; 
pave: real ; 

{ 3 

function nbh(p1 ,l,qs,qt:real):real; 
var 

x,y,u,v:real; 

begin 

x:=qs/qt; 

y:=1-2*x; 

u:=sqrt(1+y’*‘y); 

v:=u/y; 

nbh : =arctcn ( V )’*l/p1 
end; 

i — — > 

function bhdKpl ,pf ,pg,cl ,c2,c3 ,c4,ho,hb,l: real): real; 

var x,y,z:real; 

begin 

x: =si n ( (p1*ho)/(2 . 0*1) )/cos ( ( p1*ho)/(2 . 0*1 ) ) ; 
y:=ln((c4-c3)*x +sqrt(c4*c4-c3*c3)) 

-ln( (c4-c3)*x-sqrt(c4’*fe4-c3*c3) ) : 
2:5;y*((cl’*1::4-c2*c3)*l)/(c4*sqrt(c4*c4-c3*c3)*hb*pi ) ; 
bhdl : =pf- (pf-pg)*(c2/c4-2 ) 
end: 

I — — > 

procedure I:rfid2(pb,p1 ,pf,pg,c1,c2,c3,c4,ho,hb,l:real); 

var x,y, 2 :real; 

begin 

X : =s1n(p1*ho/(2. 0*1) )/cos ( pi*ho/(2 . 0*1) ) : 
y : = 1 . 0-2 . 0/pi*arctan ( (c3-c4 )*x/sqrt (c3*c3-c4*c4 ) ) ; 

2 : = (cl*c4-c2*c3 )*l/(c4*sqrt (c3’^:3-c4*c4)*hb) ; 

pb: =pf - ( pf-pQ)*(c2/c4+y*2 ) : 

end; 


begin 
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vg:- 1/pg: 

dk‘lt^i:=»33.0; 

d:-(h^-hfeed)/(hf-hfe€*d); 

{ wrlteln('U»€‘ value of d Is'.d);} 
xe1:=0.357; 

{wrIteInC ' the value of xe[k} 1s‘,xelk]);} 

ml Ck]: = (pof[k]*2.t781*civ*nob[k]*0. 138*0. 138*l)/(xetk3*hfg+deltah); 

{mg/xelkl;) 
po:=2/(v1+vf ); 
qt:=hf+(xe[k3’*hfg)-h1 ; 
qs:=deltah; 

cl :=qt/(2*hfg)- (hf-h1 )A»^gj 
c2:=-qt/(2*hfg) : 
si:=pg/pf*s; 
hyd: =0.0 1024; 

A1:=0. 01907; 
a2: =0.008827; 
g:=mUk]/nob[k3; 

{wrltelnCvalu of g is’,g);} 

gun1t:=g/a2 : 

vl:=0.95; 

vll -.=0.95*9/6.0; 

reyn : =hyd*pf*v 1/mul ; 

rey n 1 : =hyd*pf*v 1 1 /mul ; 

alphaeU=xeVtsi+( 1 . 0-s1 )*xel ) ; 

alph2e:=xe[k3/(s1+(1.0-si)*xeCk3); 

pe: =alphae*pg+ ( 1 .0-al|:rfiae)*pf ; 

pel :=alphael*pg+ ( 1 .0-alfrfiael )*pf ; 

f rf : =0 . 1 84/exp(0 . 2*ln ( reyn ) ) ; 

f rf 1:=0. 184/exp(0.2*ln(r6ynl ) ) ; 

R_: = 1 . 0/3 . 0* ( 1 . 0+ 1 . 0/ ( 1 . 0-alphae )+sqr ( 1 . 0/ ( 1 . 0-alphae) ) ) ; 
ra-;=vf*(sqr(l .0-xeCk3)/(l .0-alphae)+sqr(xe[k3)/alphae*vg/vf-l .0) ; 
accprd:=rac*sqr(gunit) ; 

sppad-. =-0 . 7* ( 1 , 0/Al*Al - 1 . 0/ (a2*a2 ) )*g/pf*g : 
tpprd:=( 1 .0/(A1*a2)- 1 .0/(Al*Al ) )*^g*(sqr( 1 . 0-xeCk3 ) 

/(pf*( 1 .0-alphae) )+xetk3*xe[k3/(p^alphae) ) ; 
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c?.: -S.H( 1 .0-si )"'ol ; 

C^:=(1-Sl )*C?; 

ho:= nb}i(p1 ,l,qs,qt); 
hb:=l-ho: 

{writelnCthe value of non_bo1l1ng hei^t Is’, ho);} 

{writelnC 'the value of boiling height is’.hb);} 

prf:=frf*pf*sqr(vl)*(ho+R_*hb)/(2.0*hyd): 

if sqr(c4)>sqr(c3) 
then begin 

pb:=bhdUpi ,pf ,pg,cl ,c2,c3,c4,ho,hb,l) ; 

end 

else begin 

bhd2(pb,pi ,pf .pg.cl ,c2,c3,c4,ho,hb,l) ; 

end; 

f 1 : =f rf 1 *sqr (rfi*4 . 0/ ( pi *pe 1 ’’‘dco^dco ))*pel/(2.0*dco); 
f 2 : =f rf I’^sqrCnf'A . 0/ (pf*pi* 1 0 . 4^ ) )*pf /( 2 . 0* 1 . 056 ) ; 
remprd: =accprd+sppad+tpprd+prf ; 
pave: = l/l’'‘(po*ho+pb*hb) ; 

hchCk] : = (9.81*pf*l-9.81*pave’''l-f2’''l-remprd)/(f l+f2-9-8l’*‘Pf'’^P®*^*®l ^ ' 
end; 

begin{ma1n program} 
xeCl]t=0.33; 
for k:=2 to ? do 
xe[k3:=0-30; 
nobtl]:=60; 
nob[Z]:=104; 
nob[3]r = n2; 
nobt4l:=176; 
nobC5]:=224: 
nob[63:=100; 
nobC?3:=96; 
poft13‘.=0.988; 
pof[23-.=0.92; 
pof L33:=0.84; 
pof [43 -.=0.7652; 
pof[53:= 0.685; 
pof[63: =0.567; 
pof C73:*0.44; 
for k:=l bo 1 do 
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rsrht(xe,pc>f ,r>ob,in1 ,tK;h): 
for k:*? to 7 do 
[k3*.=0.5; 
for k:*2 to 7 do 
begin 

hrf>ile(hch[k]<=hch[1]) do 
begin 

rsrht(xe,pof ,nob,m1 ,hcfO ; 
xe[k] :=xe[k]+0.001 
end 

end; 

incum:=0.0; 
xcLHn:’=0.0; 
k: = 1; 

wrfi11e(k<=7) do 
begin 

mcum: =mcum+m1 [k 3 ; 
xcum:=xcum+ml tk3"'xe[k3 ; 
gl [k3:=ni1 [k3/nc^Ck3 ; 

wrltelnC ' ,gUk3: 10:2, ’ * ,xe[k3-. 10:2 , ' ' ,hch[k3: 10:2) ; 

k:®k+l 

end; 

niaw:=mcum/872; 

xav : =xcLBn/(mav*872 ) ; 

wr1teln( ' valu of xav is' , xav: 15:2) : 

writelnC'valu of roav is' ,mav: 15:2): 

^d, {main program) 
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